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Abstract. Recently, it has been proposed that the natural connectivity can be used 
to efficiently characterise the robustness of complex networks. Natural connectivity 
quantifies the redundancy of alternative routes in a network by evaluating the weighted 
number of closed walks of all lengths and can be regarded as the average eigenvalue 
obtained from the graph spectrum. In this article, we explore the natural connectivity 
of random graphs both analytically and numerically and show that it increases linearly 
with the average degree. By comparing with regular ring lattices and random regular 
graphs, we show that random graphs are more robust than random regular graphs; 
however, the relationship between random graphs and regular ring lattices depends on 
the average degree and graph size. We derive the critical graph size as a function of 
the average degree, which can be predicted by our analytical results. When the graph 
size is less than the critical value, random graphs are more robust than regular ring 
lattices, whereas regular ring lattices are more robust than random graphs when the 
graph size is greater than the critical value. 
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1. Introduction 

Networks are everywhere. Many systems in nature and society can be described 
as complex networks. Examples of networks include the Internet [I], metabolic 
networks [2], electric power grids [5J, supply chains [4], urban road networks [5], world 
trade web j6] and many others. Complex networks, more generally, complex systems 
have become pervasive in today's science and technology scenario and have recently 
become one of the most popular topics within the interdisciplinary area involving 
physics, mathematics, biology, social sciences, informatics, and other theoretical and 
applied sciences(see [THS]). Complex networks rely for their function and performance 
on their robustness, that is, the ability to endure threats and survive accidental events. 
Due to their broad range of applications, the attack robustness of complex networks has 
received growing attention. 

Simple and effective measures of robustness are essential for the study of robustness. 
A variety of measures, based on different heuristics, have been proposed to quantify 
the robustness of networks. For instance, the vertex (edge) connectivity of a graph 
is an important, and probably the earliest, measure of robustness of a network [10J. 
However, the vertex (edge) connectivity only partly reflects the ability of graphs to 
retain connectedness after vertex (or edge) deletion. Other improved measures include 
super connectivity [XT] , conditional connectivity [12], restricted connectivity [15] . fault 
diameter [H], toughness [13], scattering number [IB] , tenacity [IT], the expansion 
parameter \18\ and the isoperimetric number [19]. In contrast to vertex (edge) 
connectivity, these new measures consider both the cost to damage a network and how 
badly the network is damaged. Unfortunately, from an algorithmic point of view, the 
problem of calculating these measures for general graphs is NP-complete. This implies 
that these measures are of no great practical use within the context of complex networks. 
Another remarkable measure used to unfold the robustness of a network is the second 
smallest (first non-zero) eigenvalue of the Laplacian matrix, also known as the algebraic 
connectivity. Fiedler |20j showed that the magnitude of the algebraic connectivity 
reflects how well connected the overall graph is; the larger the algebraic connectivity 
is, the more difficult it is to cut a graph into independent components. However, the 
algebraic connectivity is equal to zero for all disconnected networks. Therefore, it is too 
coarse a measure to be used for complex networks.. 

An alternative formulation of robustness within the context of complex networks 
emerged from the random graph theory [21] and was stimulated by the work of Albert et 
al. [22]. Instead of a strict extreme property, they proposed a statistical measure, that 
is, the critical removal fraction of vertices (edges) for the disintegration of a network, 
to characterise the robustness of complex networks. The disintegration of networks can 
be observed from the decrease of network performance. The most common performance 
measurements include the diameter, the size of the largest component, the average path 
length, the efficiency [23] and the number of reachable vertex pairs [24] . As the fraction 
of removed vertices (or edges) increases, the performance of the network will eventually 
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collapse at a critical fraction. Although we can obtain the analytical critical removal 
fraction for some special networks [251429] . generally, this measure can only be calculated 
through simulations. 

Recently, Wu et al. |30j showed that the concept of natural connectivity can be 
used to characterize the robustness of complex networks. The concept of natural 
connectivity is based on the Estrada index of a graph, which has been proposed 
to characterize molecular structure [31J, bipartivity [32J, subgraph centrality [33] 
and expansibility [34"|l35]. Natural connectivity has an intuitive physical meaning 
and a simple mathematical formulation. Physically, it characterises the redundancy 
of alternative paths by quantifying the weighted number of closed walks of all 
lengths leading to a measure that works in both connected and disconnected graphs. 
Mathematically, the natural connectivity is obtained from the graph spectrum as an 
average eigenvalue and it increases strictly monotonically with the addition of edges. 
Abundant information about the topology and dynamical processes can be extracted 
from a spectral analysis of the networks. Natural connectivity sets up a bridge 
between the graph spectra and the robustness of complex networks, allowing a precise 
quantitative analysis of the network robustness. 

In our previous study [36], we have shown that the natural connectivity of regular 
ring lattices is independent of the network size and increases linearly with the average 
degree. In this study, we investigate the natural connectivity of random graphs and 
compare it with regular graphs. The article is structured as follows. In Section 2, 
we introduce the concept of natural connectivity and some basic elements of random 
graphs. In Section 3, we derive the natural connectivity of random graphs. In Section 
4, we compare the natural connectivity of random graphs with that of regular graphs. 
Finally, the conclusions are presented in Section 5. 

2. Preliminaries 

2.1. Graph and Natural Connectivity 

A complex network can be viewed as a simple undirected graph G(V,E), where V is 
the set of vertices, and E C V x V is the set of edges. Let N — \V\ and M — \E\ be 
the number of vertices and the number of edges, respectively. Let A{G) = (ay)jvxJV be 
the adjacency matrix of G, where 1 if vertices Vi and Vj are adjacent, and 

= otherwise. It is obvious that A(G) is a real symmetric matrix. We thus 
let Ai > A 2 > ... > Xn denote the eigenvalues of A which are usually also called the 
eigenvalues of the graph G itself. The set {Ai, A 2 , ...Ajv} is called the spectrum of G. 
The spectral density of G is defined as the sum of the 5 functions as follows 



which converges to a continuous function when N — )■ oo, where <5(A— Aj) = 1 if A=Aj; 
and 5(X— Aj) = 0, otherwise. 




(1) 
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A walk of length A; in a graph G is an alternating sequence of vertices and edges 
v eiVie 2 .-.ekVk, where Vi £ V and ej = (t>j_i,t>j) G E. A walk is closed if t> = w fc- 
The number of closed walks is an important index for complex networks. Recently, 
we have defined the redundancy of alternative paths as the number of closed walks 
of all lengths [30]. Considering that shorter closed walks have more influence on the 
redundancy of alternative paths than longer closed walks and to avoid the number of 
closed walks of all lengths to diverge, we scale the contribution of closed walks to the 
redundancy of alternative paths by dividing them by the factorial of the length k. That 
is, we define a weighted sum of numbers of closed walks S = J2T=o n k/k\, where n k is 
the number of closed walks of length k. This scaling ensures that the weighted sum does 
not diverge and it also means that S can be obtained from the powers of the adjacency 
matrix: 

n k = trace(A k )=J2^ (2) 
i=i 

Using Eq. [2, we can obtain 

oo oo N \k N oo \k N 

A;=0 fc=0i=l i=lk=0 ^ i=\ 

Hence the proposed weighted sum of closed walks of all lengths can be derived 
from the graph spectrum. We remark that Eq. [3] corresponds to the Estrada Index of 
the graph [31], which has been used in several contexts in the graph theory, including 
bipartivity [32J and subgraph centrality [33]. The natural connectivity can be defined 
as the average eigenvalue of the graph as follows. 

Definition [30] Let A(G) be the adjacency matrix of G and let Ai > A2 > ... > Xn 
be the eigenvalues of A(G). Then the natural connectivity or natural eigenvalue of G is 
defined by 

A = ln(xy7iv) (4) 
It is evident from Eq. [4] that A at < A < Ai- 

A regular ring lattice RRLn 2k is a IK — regular graph with N vertices on a one- 
dimensional lattice, in which each vertex is connected to its 2K neighbors (K on either 
side). In a previous study |36j, we have investigated the natural connectivity of regular 
ring lattices and shown that random regular graphs are less robust than regular ring 
lattices based on natural connectivity. 

Theorem 2.1. J3b^ Let RRL^^k be a regular ring lattice. Then the natural connectivity 
of RRL N 2 k is 

Xr Ni2K = In (l (2XZ2) I +o(l) (5) 
where o(l) — > as N — > 00. 
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2.2. Erdos-Renyi Random Graphs 

Random graphs have long been used for modelling the topology of systems made up of 
large assemblies of similar units. The theory of random graphs was introduced by Erdos 
et al. [37]. A detailed review of random graphs can be found in the classic book |21j. 
A random graph is obtained by starting with a set of iV vertices and adding edges 
between them at random. In this article, we study the random graphs of the classic 
Erdos-Renyi model Gn, p , in which each of the possible C% = N(N — l)/2 edges occurs 
independently with probability p. Consequently, the total number of edges M is a 
random variable with the expectation value E(M) = p-C 2 N and then the average degree 
< k >= (N - l)p « Np. 

Random graph theory studies the properties of the probability space associated 
with graphs with iV vertices as iV — > oo. Many properties of such random graphs can 
be determined using probabilistic arguments. We say that a graph property Q holds 
almost surely for Gn, p if the probability that Gn, p has property Q tends to one as 
iV — y oo. Furthermore, Erdos et al. described the behavior of Gn, p very precisely for 
various values of p [38J. Their results showed that: 

(i) If Np < 1, then a graph G^, p will almost surely have no connected components of 
size larger than o(hxN); If Np > 1, then a graph Gn, p will almost surely have a 
unique "giant" component containing a positive fraction of the vertices. 

(ii) If Np < InN, then a graph Gn, p will almost surely not be connected; If Np > IniV, 
then a graph Gtv iP will almost surely be connected. 

It is well known that the largest eigenvalue Ai of Gn, p is almost surely Np[l + o(l)] 
provided that Np ^> lniV(see [391110]). Moreover, according to the famous Wigner's 
law or semicircle law [41 J, as N — > oo, the spectral density of Gn, p converges to a 
semicircular distribution as follows 



P(A) = 







< R 


I o 


A 


> R 



(6) 



where R = 2 J Np(l — p) is the radius of the "bulk" part of the spectrum. 



3. Natural connectivity of random graphs 

When iV — > oo, with continuous approximation for Aj, Eq. [4J can be rewritten in the 
following spectral density form 

A = In p(\)e x d\^ = In (M A (1)) (7) 

where p(A) is the spectral density and M\(t) is the moment generating function of p(A). 
Consequently, we obtain the natural connectivity of Erdos-Renyi random graphs with 
p > In N/N 

A = In ( ( +R p(\)e x d\ + e Al /iV] = In (m a (1) + e Np /N) (8) 
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-R2y/R 2 -\ 2 Aj , 2 r+RVR 2 -X 2 x 



Substituting A = Rcos(9) into Eq. (8), we obtain that 

M A (1) = - e Rcos{e) sin 2 (9)d9 (10) 



7T JO 

Note that 



I a ( x ) = ^ /2) " / " e xcos{e) sin 2a (9)d9 (11) 



7f 1 / 2 T(a + 1/2) Jo 

where 7 a (x) is the modified Bessel function and T(x) is the Gamma function. Then we 
obtain that 

I,(R) = - T e Rcos{e) sin 2 (9)d9 (12) 



7T JO 

Using Eq. (11), we can simplify Eq. (9) as 

M x (l)=2h(R)/R (13) 
Substituting Eq. (12) into Eq. (7), we obtain that 

^K^+^H- inW +K i+ ^) (i4) 

Now we propose two lemmas first. 

Lemma 3.1. As N — > oo, f(p) = 2NIi(R)/(e Np R) is a monotonically decreasing 
function for In N/N < p < 1 — In N/N, where R = 2J Np(l — p) . 



Proof. It is easy to know that 2y / lniV(l - In N/N) < R < VN for In N/N < p < 
1 — In N/N. Then as N — > oo, we have R — > oo. Note that, for the large values 
of x 3> \a 2 — 1/4 1, the modified Bessel functions I a (x) have the following asymptotic 
forms [32] 

IJx) -> -^=e x (15) 
V 2irx 

Thus, for In N/N < p < 1 - In JV/iV, we obtain 

AW "T^e* (16) 



Then we have 



Note that, 



fW^^I--^ ( 17 ) 



rffb) > e^-^(^-iv)/;3/ 2 -| fi i/ 2 .^. e H-^ 

dp ~^ R 3 



= f V f (iV(2 - 4 P - i?) - ^iizM) < o 

Therefore, we prove that, as iV — > oo, /(p) is a monotonically decreasing function for 
\nN/N < p < 1 - \nN/N. □ 
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Lemma 3.2. Let p c = N e 1 lnA r , where < £ < 1. Then we have f(p c ) — > as 
N -»■ oo. 

Proof. Note that 0<£<Cl, thus we have p c — )■ and 1— p c — ^ 1 as iV - >■ oo. Then we 
obtain that R Pc — > 2a/ N £ In N. Therefore, we have 



fl£' c ' V t ( 2 VAf £ In AT) 

Afl-3e/4 



20F ' (in AT) 3/4 2v^F (In AT) 3 / 4 



(19) 



jyl-3e/4 e -(ViV e InJV-l) +1 

= 2v^F ' (i n at)3/4 

Since -\/iV £ In iV > 1 as -> oo, we obtain that 

jyl-3e/4 -( VlVMnlV-l) 2 + l e jyl-3e/4-Af e 
/(p c ) -»■ — -= TT A PS tjj -> (20) 

V ; 2v^F (lniV) 3/4 2^/?(lniV) 3/4 V ^ 

The proof is completed. □ 

From Lemmas 3.1 and 3.2, it is easy to derive that, for p c < p < 1 — p c , 
f(p) < f(Pc) — t as — > oo. Consequently, we obtain the following theorem. 

Theorem 3.3. Let Gn jP be a random graph with N e ~ l \nN < p < 1 — N 6 ^ 1 In N , where 
0<£<1. Then the natural connectivity of Gjv iP almost surely is 

X = Np — ln(iV) + o(l) (21) 

where o(l) -} as iV 4 oo. 

From Eq. (20), we know that the natural connectivity of random graphs increases 
linearly with edge density p given the graph size N. Note that < k >= Np; thus, 
we also observe that the natural connectivity of random graphs increases linearly with 
the average degree given the graph size N. To verify our result, we simulate 1000 
independent Gn. p and compute the average natural connectivity for each combination 
of TV and p. In Fig. 1, we plot the natural connectivity of random graphs with both 
numerical results and analytical results. We observe that the numerical results agree 
well with the analytical results. 



4. Comparisons between random graphs and regular graphs 

We generate Erdos-Renyi random graphs Gn, p , regular ring lattices RRL^,2K, and 
random regular graphs RRGn^k using the algorithm in [H]. We compare the natural 
connectivity of random graphs with two other types of regular graphs with the same 
number of vertices and edges, i.e., p ps 2K/N and < k >= 2K. The results are shown 
in Fig . 2. We find that regular ring lattices and random graphs are always robustness 
than random regular graphs. However, the curves of regular ring lattices cross those of 
random graphs; furthermore, random graphs are more robust than regular ring lattices 
prior to crossings (dense networks), whereas regular ring lattices are more robust than 




Figure 1. Natural connectivity of random graphs: (a) A vs. p with N = 100 (circles), 
500 (triangles) and 1000 (diamonds); (b) A vs. N with p — 0.1 (circles), 0.3 (triangles), 
0.5 (diamonds). Each quantity is an average over 1000 realizations. The lines represent 
the corresponding analytical results according to Eq. (20). 




random graphs over crossings (sparse networks). This means that there is a critical 
graph size N c , that is as a function of K. For example, for K = 5, we find that iV c 60. 

For large values of K, we can analytically predict the values of N c using Eq. (4) 
and Eq. (20) as follows 

k \ K 

= Np-ki(N) = 2K-\n(N) => N c « e 2 ^- /o(2 ' 2 ' - 2) (22) 

The results are shown in Fig. 3. Moreover, we also find that there is a critical value p c 
or K c as a function of graph size N. Regular ring lattices are more robust than random 
graphs when the edge density p < p c , whereas random graphs are more robust than 
regular ring lattices when the edge density p > p c . 

To explore the critical behaviors of graphs in depth, we randomise regular ring 
lattices by random rewiring [15] and by random degree-preserving rewiring [16], which 
leads to random graphs and random regular graphs, respectively. In Fig. 4, the natural 
connectivity is represented as a function of the number of rewirings, starting from regular 
ring lattices with iV = 30 < N c and iV = 100 > N c , where K = 5. We find that the 
natural connectivity decreases during the process of random degree-preserving rewiring 
and equals to the value of a random regular graph finally. It means that regular ring 
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Figure 2. Natural connectivity of random graphs of regular ring lattices RRLn, 2 k 
(squares), random regular graphs RRGn,2K (triangles) and random graphs Gn, p 
(diamonds) with the same number of vertices and edges. From bottom to top, the 
symbols correspond to K = 3, 4, 5, respectively. Each quantity is an average over 1000 
realizations. 
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Figure 3. The critical value N c as a function of graph size K according to Eq. (21). 
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lattices are more robust than random regular graphs for both N = 30 and iV = 100. The 
case of random rewiring is more complicated. Different processes of random rewiring 
for N = 30 and iV = 100 are shown in Fig. 4. The natural connectivity increases 
during the process of random rewiring for N = 30 < N c ; however, for N = 100 > N c , 
the natural connectivity first decreases during the process of random rewiring and then 
increases during the process of random rewiring; finally, equals to the value of a random 
graph finally (smaller than the value of a regular ring lattice). It means that randomness 
increases the robustness of a dense regular ring lattice, but decreases the robustness of 
a sparse regular ring lattice. 
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Figure 4. Natural connectivity during the processes of random rewiring (diamonds) 
and random degree-preserving rewiring (triangles) starting from regular ring lattices 
with N — 30 (a), N = 100 (b), where K — 5. The solid lines represent the values of 
random regular graphs and the dashed lines represent the values of random graphs. 
Each quantity is an average over 1000 realizations. 



5. Conclusions 

We have investigated the natural connectivity of Erdos-Renyi random graphs Gjy, P in this 
article. We have presented the spectral density form of natural connectivity and derived 
the natural connectivity of random graphs analytically using the Wigner's semicircle law. 
In addition, we have shown that the natural connectivity of random graphs increases 
linearly with edge density p given a large graph size N. The analytical results agree 
with the numerical results very well. 
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We have compared the natural connectivity of random graphs Gat iP with regular 
ring lattices RRL^ 2 k and random regular graphs RRGn^k with the same number of 
vertices and edges. We have shown that random graphs are more robust than random 
regular graphs; however the relationship between random graphs and regular ring lattices 
depends on the graph size N and the edge density p or the average degree < k >. We 
have observed that the critical value N c as a function of K, and the critical value p c 
and K c as a function of graph size N, which can be predicted by our analytical results. 
We have explored the critical behavior by random rewiring from regular ring lattices. 
We have shown that randomness increases the robustness of a dense regular ring lattice, 
but decreases the robustness of a sparse regular ring lattice. Our results will be of great 
theoretical and practical significance to the network robustness design and optimization. 
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